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We study a supersymmetric model for strongly interacting lattice fermions in the presence of 
a staggering parameter. The staggering is introduced as a tunable parameter in the manifestly 
supersymmetric Hamiltonian. We obtain analytic expressions for the ground states in the limit of 
small and large staggering for the model on the class of doubly decorated lattices. On this type of 
lattice there arc two ground states, each with a different density. In one limit wo find these ground 
states to be a simple Wigner crystal and a valence bond solid (VBS) state. In the other limit we 
find two types of quantum liquids. As a special case, we investigate the quantum liquid state on the 
one dimensional chain in detail. It is characterized by a massless kink that separates two types of 
order. 
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I. INTRODUCTION 

A central theme in modern day solid state physics is the study of the electronic properties of materials in which the 
electrons are free to move, but subject to (strong) repulsive interactions. In this regime one typically faces the difficulty 
that perturbative tools are unreliable. Consequently, even at the level of relatively simple model Hamiltonians, the 
behavior of such systems is notoriously difficult to analyze. This challenge motivates the study of a supersymmetric 
model for strongly interacting, spinless fermions [T]. In these models a judicious tuning of kinetic and potential terms 
leads to systems which possess supersymmetry. These models show various interesting features that are closely linked 
to the supersymmetry. First, supersymmetry provides a rich mathematical structure that allows for a considerable 
degree of analytic control in the regime where more standard perturbative techniques would fail. Second, it induces 
a delicate balance between kinetic and potential terms, which results in a strong form of quantum charge frustration. 
This so-called 'superfrustration', which typically occurs for two (or higher) dimensional lattices, is characterized by 
an extensive ground state entropy 1^ . 

The degrees of freedom in the supersymmetric model are itinerant hard-core spin-less fermions. For spin-less 
fermions double occupancy is forbidden by the Pauli exclusion principle. The hard-core constraint which imposes 
nearest neighbor exclusion extends this repulsive character. At low enough densities a system of such particles 
exhibits a metallic phase described by spin-less Fermi liquid theory. On the square lattice the phase at half-filling 
(one particle every two sites) is trivial; there are two ordered patterns in which the fermions are completely jammed. 
However, for small hole-doping there is strong evidence that the hard-core spin-less fermions exhibit a stripe phase, 
where domain walls between the different orders favor being parallel [JHB]. 

It is clear that the intermediate density regime of this strongly-interacting system is non-perturbative; one has 
to account for both the kinetic as well as the potential terms. It turns out that the supersymmetric model, which 
incorporates a next-nearest neighbor repulsion for the itinerant hard-core spin-less fermions, provides a way to study 
precisely this regime. Exact results for this model on a variety of two (and higher) dimensional lattices have been 
obtained [2, J|M1] (for reviews see [HITl]). 

A key feature of the model is superfrustration [2]: a ground state degeneracy that is exponential in the number 
of lattice sites. One often finds that these ground states occur in a whole region of particle densities. For example, 
for the supersymmetric model on the square, triangular and honeycomb lattice there exist zero energy ground states 
for all rational densities v within the range [1/5,1/4], [1/7, 1/5] and [1/4, 5/18], respectively l9, . Note that at these 
densities the model is deep inside the non-perturbative regime. 

A qualitative understanding of superfrustration stems from the 3-rule, which says that particles prefer to be mostly 
3 sites apart [2] . In general it turns out that there is an exponential number of ways to satisfy the 3-rule. For certain 
lattices, however, the frustration can be dramatically reduced. These lattices have a structure that nicely accomodates 
the 3-rule. 

On the square lattice the frustration leads to a sub-extensive ground state entropy; the ground state degeneracy is 
exponential in the linear size of the system. This result is for the square lattice wrapped around the torus, for the 
cylinder and the plane there is compelling evidence that the system supports quantum critical edge modes |10| . 

Quantum criticality is also observed for the supersymmetric model in one dimension, where the system turns out 
to be integrable. The continuum theory is a Luttinger liquid with special values of the Luttinger parameters, such 
that the theory is supersymmetric [3 [TS] . In recent work [ini [TT] , it was shown that the one dimensional model also 
enjoys special features, such as scale-free properties, which are again directly related to supersymmetry. 

In this paper we study the supersymmetric model extended with a tunable parameter, which for reasons that will 
become clear later, we call the staggering parameter pi lTSHTS] . We show that for the lattices with reduced frustration 
we can obtain analytic expressions for the zero energy states in the limits where this parameter is very large and very 
small. 

Interestingly, we find a whole class of graphs where, in the limit of a very small staggering parameter, the ground 
state is the equal weight superposition of all states in the low energy sector. Such liquid states typically arise in the 
context of quantum dimer models, where this is a consequence of the Perron-Frobenius theorem. Famous models with 
liquid ground states are Kitaev's toric code [TH] and Levin and Wen's string-nets pUj, which exhibit topological order, 
and the Rokhsar-Kivelson (RK) model ^T], which exhibits criticality. The class of models we discuss in this paper, 
however, seem to be of a different type. We find that for a particular choice of lattice, our effective Hamiltonian 
precisely maps to the constrained quantum-clock model that arises in the context of bosons in a tilted optical lattice 
[22]. For this model it has been shown that the ground state has no long range order, which suggests that it is in a 
gapped phase. In our model this quantum liquid state has one particle per plaquette. We find that there is another 
quantum liquid ground state at a higher density, which has one particle per link. The two quantum liquid states arise 
from a projection of a product state of single particle states with one particle per vertex plaquette and one particle 
per link respectively. The procedure is reminiscent of the construction of the Affleck-Kennedy-Lieb-Tasaki (AKLT) 
state [221 121] ■ Besides the quantum liquids, we will find that the model also supports other gapped phases. One being 
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a charge ordered crystalline phase, the other a valence bond solid (VBS) that does not break any lattice symmetry. 
In contrast to the quantum liquids, the ground states in the latter two phases are simple product states. 

As a special case of the class of models with liquid ground states we discuss the staggered supersymmetric model 
on the one dimensional chain with open boundary conditions. It is well known [TlllSj that the supersymmetric model 
on the chain without staggering is quantum critical. In [16j it was shown that the staggering parameter corresponds 
to a relevant operator in the continuum theory. In this paper we show that in one limit of the staggering parameter, 
the ground state of the chain with open boundary conditions is a liquid state. This liquid state contains a delocalized 
massless kink. We will see that, as a consequence, the spectrum in that limit has a zero dimensional massless character. 

In the last part of the paper, we speculate that for lattices that do not have reduced frustration an analysis in the 
same spirit as the one presented here may lead to the discovery of other exotic phases. There are indications for RK 
type phases and the existence of one dimensional critical modes in two dimensional lattices. 

In the remainder of this introduction we first discuss the basics of supersymmetry for lattice fermions. We then 
define the model and discuss the consequences of the staggering parameter. Finally, we introduce cohomology. This 
mathematical technique has proven extremely powerful in the analysis of the supersymmetric model and will also be 
essential for the analyses in this work. 

A. Supersymmetry for lattice fermions 

Supersymmetry is a symmetry between fermionic and bosonic degrees of freedom (see [25j for a general reference) . 
It plays an important role in theoretical high energy physics, where various theories that go beyond the standard 
model require supersymmetry for a consistent formulation. In these theories all the known elementary particles are 
accompanied by yet to be discovered superpartners. 

In the lattice model discussed here the physical particles are spinless lattice fermions and the supersymmetry 
relates fermionic and bosonic many particle states with an odd and even number of the lattice fermions respectively. 
In the Af — 2 supersymmetric theories that we consider, a central role is played by the operators Q and , called 
supercharges, which have the following properties [261 

• Q adds one fermion to the system and takes out one fermion from the system. 

• The supercharges are fermionic operators and thus nilpotent: = (Q^)^ = 0. 

• The hamiltonian is the anti-commutator of the supercharges, H — {Q, Q^}, and as a consequence it commutes 
with the supercharges and conserves the number of fermions in the system. 

Imposing this structure has some immediate consequences; supersymmetric theories are characterized by a positive 
definite energy spectrum and a twofold degeneracy of each non-zero energy level. The two states with the same energy 
are called superpartners and are related by the supercharge. 

From its definition it follows directly that H is positive definite: 

= iQiV')i' + igV>i'>o. 

Furthermore, both Q and commute with the Hamiltonian, which gives rise to the twofold degeneracy in the energy 
spectrum. In other words, all eigenstates with an energy Eg > form doublet representations of the supersymmetry 
algebra. A doublet consists of two states \s),Q\s), such that Q^s) = 0. Finally, all states with zero energy must be 
singlets: Q\g) — Q^g) = and conversely, all singlets must be zero energy states [22]. In addition to supersymmetry 
our models also have a fcrmion-number symmetry generated by the operator F with 

[F, gt] = -Qt and [F, Q] = Q. 

Consequently, F commutes with the Hamiltonian. From the commutators of F with the supercharges, one finds that 
Q and change the number of fermions by plus or minus one unit, so that the supercharges map bosonic states to 
fermionic states and vice versa. 

B. The model 

We now make things concrete and define supersymmetric models for spin-less fermions on a lattice or graph with L 
sites in any dimension, following [T. The operator that creates a fermion on site i is written as cl with {cj, cj} = 6ij. 
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To obtain a non-trivial Hamiltonian, we dress the fermion with a projection operator: P^,;^ = Ylj next to i(l ~ '^j'^j): 
which requires all sites adjacent to site i to be empty. We now define the supercharges as 

Q = 5]A*4P(,), =Y,^'^C,P{^)■ (1) 

One can check that the supercharges are nilpotent for any complex numbers A^. It follows that the Hamiltonian of 
these hard-core fermions reads 

i j next to i i 

The second term contains a next-nearest neighbor repulsion, a chemical potential and a constant. The details of the 
latter terms will depend on the lattice we choose. The first term is just a nearest neighbor hopping term for hard-core 
fermions. The projection operators in this term ensure that the nearest neighbor exclusion is obeyed by the hard-core 
fermions. Another way to impose this condition would be to include a nearest neighbor repulsion of strength V and 
send V to infinity. 

Physical systems which realize this type of interactions are ultra cold gasses of Rydberg atoms. The Van der Waals 
interaction between Rydberg atoms falls off as l/r^, where r is the distance between two Ryberg atoms. It follows 
that the nearest neighbor repulsion is almost two orders of magnitude larger than the next-nearest neighbor repulsion. 
This results in what is called a Rydberg excitation blockade; atoms within a certain distance of a Rydberg atom, 
cannot be excited to the Rydberg state. It follows that Rydberg atoms can be accurately modelled by hard-core 
fermions. 

In this paper we focus on lattices S which naturally break up into two sublattices Si and S2 = S\Si. On these 
lattices we choose the following staggering: 

A. = |M"'^? (2) 
[ y for I e D2 

It follows that the hopping amplitudes between two 5*1 sites, an 5*1 and an 5*2 site and two S2 sites are 1, y and 
respectively. Furthermore, we find that the chemical potential on the sublattice Si is independent of y, whereas the 
chemical potential on 5*2 is proportional to y^. The consequences of the staggering on the interaction terms is a bit 
more subtle and will depend on the details of the lattice. 



C. Cohomology 

For the supersymmetric models, cohomology has proven to be a very powerful tool to extract information about 
the zero energy ground state(s) of the models (see for example [U El EHUl E] and references therein). The key 
ingredient is the fact that ground states are singlets; they are annihilated both by Q and . This means that a 
ground state \g) is in the kernel of Q: Q\g) — and not in the image of Q, because if we could write \g) = Q\f), then 
(I/), \g)), would be a doublet. So the ground states span a subspace Hq of the Hilbert space H of states, such that 
Hq — kerQ/Im Q. This is precisely the definition of the cohomology of Q. So the ground states of a supersymmetric 
theory are in one-to-one correspondence with the cohomology of Q. It follows that the solution of the cohomology 
problem gives the number of zero energy states for each particle number sector. 

We compute the cohomology using the 'tic-tac-toe' lemma of This says that under certain conditions, the 

cohomology Hq for Q = Qi + Q2 is the same as the cohomology of Qi acting on the cohomology of Q2- In an 
equation, Hq = 'Hq^(T-Lq^) = 'H12, where Qi and Q2 act on different sublattices 5*1 and S'2. We find 'H12 by first 
fixing the configuration on all sites of the sublattice 5*1, and computing the cohomology T-Lq^- Then one computes the 
cohomology of Qi, acting not on the full space of states, but only on the classes in T-Lq^. A sufficient condition for 
the lemma to hold is that all non-trivial elements of H12 have the same /2 (the fermion- number on S'2). 

Although ground states are in one-to-one correspondence with cohomology elements, the two are not equal unless 
the cohomology element happens to be a harmonic representation of the cohomology. Harmonic representations are 
elements of the cohomology of both Q and . So they are annihilated by both supercharges, which is precisely the 
property of a zero energy state. It follows that, although the solution of the cohomology problem gives the number 
of ground states, it typically does not give the ground states themselves. In this paper, however, we show that for 
the lattices with reduced frustration the cohomology computation gives the exact ground states in the limits where 
the staggering is very large and very small. This property was also successfully exploited in [16l [TT] . The fact that 
cohomology elements can provide more information about the ground states than just their degeneracy and particle 
density, was also observed in [101 HH HI] • 
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II. COHOMOLOGY FOR STAGGERED LATTICES 



In this paper we investigate the properties of the supersymmetric model in the Umits of very small and very large 
staggering. It follows that wo can treat the problem perturbatively in either y or 1/y. As a first step we compute the 
ground states to zeroth order in the two limits. When we find that the number of zeroth order ground states is larger 
than the number of ground states at finite y, the second step is to do degenerate perturbation theory. We find that 
the degeneracy is lifted at second order in perturbation theory. In both steps cohomology plays an essential role. 

First of all, we will find the zeroth order groimd states by solving the cohomology problem in the two limits. 
Remember that we introduced the staggering parameter y. such that the supercharge takes the form Q = Qi + yQ2 
and the Hamiltonian reads H = {Q,Q^} = Hq + yHi + j/^i?2. It follows that for y = we have Q = Qi and for 
y — >■ 00 we have Q « yQ2. For the lattices we consider, we can easily compute both T-Iq^ and "Hqj and we will find 
that the elements of Hq^ and Hq^ respectively are precisely the zeroth order ground states in the small and large 
staggering limit. 

Second of all, cohomology tells us that if the number of ground states in either of the two limits is larger than the 

number of ground states for finite y, this degeneracy is lifted at second order in perturbation theory. Typically, we 
will find that T-Lq^ contains many more elements than T-Lq, that is, Di > D, with Di = (Hui'Hq^ and D = dim'HQ. 
It follows that for y strictly zero there are Di zero energy states. We will now prove that Di — D of these states will 
acquire an energy of order y"^ when y is no longer strictly zero, but still small. An important condition for this to be 
true is that we have that Hq — HQ^iHg^), that is, the 'tic-tac-toe' lemma should hold. 

Let Ix) be an eigenstate of the Hamiltonian, H, with energy E{y) > 0. Furthermore, we assume that Hq\xo) — 0, 
where we defined |xo) = lim^^o |x)- It follows that |xo) G T^Qi and |xo) ^ "Hg. We will assume that Hq — Hq^CHq^). 
Since |xo) ^ "^Q-iC^Qi) '^^ write Q2IX0) = |0) and Q2IX0) = W), where |(/)) and are not both zero. 

It then follows that the energy of |xo) is 



where > 0. This proves our statement above. 

In this paper we will exploit the fact that for any finite y the degeneracy at j/ = will be lifted at second order in 
perturbation theory. As we will see the effective Hamiltonian that couples the cohomology elements of Hq^ will take 
on some simple form and the zero energy states can be found exactly. 



In this section we consider the staggered one dimensional chain with open boundary conditions and length L = 3j. 
The methods explained for this special case will carry over to a variety of two dimension lattices in the next section. 
We discuss the cohomology problem in general and in the two limits of large and small staggering. We find that for 
all finite values of the staggering parameter the ground state is unique and has particle number f = j. For y strictly 
zero we find that there are 2j + 1 zero energy states. For finite staggering this degeneracy is lifted. The effective 
Hamiltonian that couples these states is obtained in second order degenerate perturbation theory. The ground state 
of the effective Hamiltonian is an equal superposition liquid state. 

We will use the following conventions. For a chain of length L = 3j we number the sites as i = 0, . . . , 3j — 1 and 
for the staggering we choose: 



E{y) « (xol^^lxo) 

= (Xol-f^olxo) +2/(xo|-ffi|xo) +2/^(xo|-f?2|xo> 

= yixoHQiQl + Q\Q2)\xo) + y^{xo\{Q2Ql + QlQ2)\xa) 

= yiQ\\xo))^\4>')+y{Qi\xo))^4>) + y\Ql\xo))^<P')+y\Q2\xo))^cl^) 

=yH{'P'm + {m) 

= Ny\ 



III. EXACT GROUND STATES FOR STAGGERED CHAIN 




1 for i = 1 mod 3 
y otherwise. 



(3) 



A. Cohomology for staggered chain 



To compute the cohomology for the chain, one uses the fact that the cohomology of a single site that can be both 
empty and occupied, is trivial. This is equivalent to the statement that the single site chain has no zero energy states. 



6 



indeed the empty and occupied state form a doublet of energy one. For the chain of length L = 3j with open boundary 
conditions, we define the sublattice as all sites i = 1 mod 3 (remember that we label the sites as < i < L — 1). 
Sublattice S2 consists of the remaining sites. Solving Hq^ is easy: it is trivial unless the sites and L — 1 are blocked 
by occupying the sites 1 and L — 2, but now 'Hq2 is still trivial unless the sites 3 and L — 4 are blocked, etc. Continuing 
this reasoning we find there is only one non-trivial element, which has all Si sites occupied. Since there is only one 
element in Hq^, finding H12 is trivial: ^12 — "Hqj. For the same reason the 'tic-tac-toe' lemma holds and we find 
that Hq consists of one non-trivial element. It follows that there is one zero energy state with j fermions. 

The above argument holds for all values of the staggering parameter, except when y = 0. In that case Q2 drops out 
alltogether so it does not make sense to compute Hq^- Instead we have Hq — Hq-^, which we can compute directly. 
Since 5*1 is a collection of isolated sites, the cohomology of Qi vanishes, unless all 5*1 sites are blocked. One easily 
verifies that there are j -I- 1 configurations with j fermions and j configurations with j + I fermions for which all Si 
sites are blocked. If we denote the sites 3i + 1 by a square and the other sites by a circle and color occupied sites 
black, we can depict these configurations as 

• □ o •□ o • • • • □ o •□o •□ o •□ o • • • • □ o •□• 

• □ o •□ o • • • • □ o oD. •□ o •□ o oD. 

: : (4) 

• □ O oD • • • • O □ • oD* •□ • oD • • • • O □ • oD* 
oD • oD • • • • o □ • oD* 

It follows that for y strictly zero, there are 2j + 1 zero energy states. In fact, in this limit the zero energy states 
are precisely the states depicted above. For y — >• 00 the zero energy state is the state with all sites occupied: 

j-i 

|V'Gs(2/^oo)) = n4+il0>- (5) 

fc=0 

Finally, we note that the Witten index, which is the trace over the ground state manifold of (—1)^, remains 
unchanged for aU values of y: W =^ (-1)^ for y > and W = {j + + j{-iy^^ = {-^Y for y ^0- 



B. Effective Hamiltonian for finite staggering 



In this section we discuss degenerate perturbation theory in the limit of small staggering parameter y. Since the 
supercharges are linear in y the Hamiltonian takes the form H = Hq + yHi + y^H2. Here Hq is purely diagonal. 
Hi is purely off-diagonal, it is the hopping between staggered and non-staggered sites, and finally, H2 consist of a 
diagonal and an off-diagonal part, where the off-diagonal part represents hopping between the sites 3i and 3i — 1 for 
1 < i < L/3- 1. 

We know that at finite staggering there is a unique ground state with zero energy and f — j- So let us write for 
this zero energy state 

|*) = |*o)+t/|*i)+y'|*2) + ... (6) 

where |^o) ~ Si=o ^fclV'fc)- The configurations \ipk) are the cohomology elements of Hq^ with j particles as shown 
m 0. They can be defined as: 

k=0 k=i 

for < i < j. Similarly, we write for the cohomology elements of Hq^ with j + 1 particles: 

i'^')=ri4n4+2i0)> (8) 

with 0<i<j-l. 

To determine the coefficients we have to go to second order perturbation theory. To first order we have 



Ho|*i> = {Ei~Hi)\^o)- 



(9) 
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Taking the inner product on both sides with \'ipi), we find 



(10) 



fe=0 



Here we used ("0(1^1 IV'fc) = 0, since Hi is purely ofF-diagonal and only has non-zero elements between configurations 
which differ in their number of particles on the sites 3i + 1 by one. Clearly, we cannot have that all ai are zero, so we 
find El =0 and the degeneracy is not lifted. 

To find the wave function to first order we write \^i) = X^m'^mlXm), where TJolXm) = Em\Xm) with £',„ > 0. We 
then take the inner product on both sides in (|9]) with \xi)- This way we find 



E 



Ofe- 



E, 



Now to second order we have 

= {El - + {E2 ~ ff2)|*o). 

Again taking the inner product with |^/) and using = 0, we find 

= {-^AiEi - Hi)\^i) + {'^i\{E2 - i?2)|*o> 
{MHi\Xm){Xm\Hi\iJk) 



(11) 



(12) 



E 

fc=0 



E 



a-k- 



E„ 



ak{4'i\H2\'4'k) 



aiE2 



(13) 



These homogeneous algebraic equations can only be solved for the coefficients ai if the determinant of their coefficients 
vanishes. The resulting secular equation reads 



where 



\Wik - SikE2\ ^ 0, 



(14) 



(15) 



The matrix W is interpreted as an effective Hamiltonian that couples the states \ipk)- It can easily be shown that W 
takes on the simple form 



W 







1 



v 



-1 

-1 2 








-1 J I 



I 1 / 



(J + l 



i + 2 





J +2 



J + l/ 



(16) 



Note that the rows of this effective Hamiltonian all sum to zero. Consequently, the zero energy ground state is an 
equal superposition of all the IV'fc): 



i*o) = Ei0fc). 



(17) 



fc=0 



Remember that we showed in section |TTj that at second order in perturbation theory Di — D states will acquire an 
energy of order y^. In the sector with j particles we have Di = dimT^Q^ —j + l and D = dimHg = 1. From 
this it follows that the (j + 1) x {j + 1) matrix W has a unique zero energy state. Another way to see that |^'o) is 
the unique ground state is via the Perron-Frobenius theorem. Note that all off-diagonal elements of W are negative 
and that repeated application of this effective Hamiltonian connects all configurations. From these two properties it 



FIG. 1: Decorated square lattice. The sublattice 5*1 is the original square lattice, the sites are depicted as squares. 
The sublattice S2 are the added sites, indicated by circles, on the links of the original lattice. For doubly periodic 
boundary conditions the right (bottom) sites are identified with the left (top) sites. 

follows by the Perron-Frobenius theorem that the ground state is unique and nodeless. For this very simple matrix, 
using the Perron-Frobenius theorem is somewhat pedantic. In fact, as we shall discuss in section [V| we can obtain 
all eigenvalues and eigenstates of this system. However, for the two dimensional lattices that we discuss in the next 
paragraph we will see that many of the effective Hamiltonians share this property. 

We have found that the degeneracy of the states with j particles is lifted in second order perturbation theory. Due 
to supersymmetry it is now directly clear that the degeneracy of the states \(f)i) with / = j + 1 is also lifted. This is 
because supersymmetry dictates that the j states with energies y^i?2 > have superpartners with equal energy and 
/ = j ± 1. Clearly, these superpartners must be linear superpositions of the states \(t>i)^ because all other states have 
energies of order one. Since there are j such states it follows that there are no zero energy states at / = j + 1. We 
conclude that there is a unique zero energy state, which has f — j particles. 

IV. EXACT GROUND STATES FOR STAGGERED DECORATED LATTICES 

We now turn to two (or higher) dimensional lattices. In this section we apply the same approach that we used in 
the previous section for the chain to a whole class of graphs (a graph is just a lattice, but the unit cell may contain 
all sites). The graphs we consider are defined as follows. We start with an original graph A and construct the graph 
A„ from the original graph by adding n — 1 additional vertices on every link. In the following the only restriction 
on the graph A is that it does not contain disconnected subgraphs. It is clear that these decorated graphs, A„, 
naturally break up into two sublattices. We have seen that the 'tic-tac-toe' lemma is then especially powerful and in 
the following we will see how this property can be exploited. 

We will first discuss the cohomology problem for the decorated graphs in general [21 HU . We will then restrict 
our attention to the doubly decorated graphs, A3. For these graphs we will show that upon staggering we can find 
quantum liquid, crystalline and VBS type states. 

A. Cohomology for decorated graphs 

The cohomology problem for the graphs of type A3,„ can be solved for general m. For graphs of type Aam-ti there is 
no general solution to the cohomology problem, however, one can relate it to the cohomology problem of the original 
lattice A. These results were obtained on a homotopy level using Alexander dualities in ^28,, but can also be derived 
using the 'tic-tac-toe' lemma 

Let us resume the cohomology problem of the graph A3, which was discussed in [2]. In the following the original 
graph A can have open or closed boundary conditions, but we assume that no changes are made to the boundary 
conditions once the additional sites are added to the links of A. We postpone the discussion on how the results 
are modified when the boundary conditions are changed by cutting links in the A3 graph to the end of this section. 
Consider the original graph A as the subgraph Si and the additional sites as the subgraph 82- The graph S2 is a 
collection of two site chains. Remember that an isolated site that can be both empty and occupied has a trivial 
cohomology. It follows that the non-trivial elements in Hq^ have the Si sites neighboring a two site chain on 5*2 either 
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D) 
\R) 

L) 







FIG. 2: Mapping of elements of the cohomology of Qi aX F = N\ to projected clock model states. On the left we 
show the four possible states of each vertex. On the right we show, as an example, a cohomology element of Hq^ 
with F = Na for the doubly decorated square lattice with open boundary conditions. 



both empty or both occupied. It follows that there are only two non-trivial elements in Hq^, one with completely 
empty and one with Si completely filled. If we leave Si completely empty, we obtain one non-trivial element of 'Hq2 
in the sector with L\ fermions, where L\ denotes the number of links in the original graph A. The element with Si 
completely filled, clearly has fermion number N\, where A^a denotes the number of vertices in the original graph A. 
It now quickly follows that these two elements are also in ?^i2, since within Hq^ both states are in the kernel of Qi 
and not in the image of Qi. Now remember that for the 'tic-tac-toe' lemma to hold a sufficient condition was that all 
elements in ■H12 have the same number of fermions on sublattice S'2. Here this condition is clearly not met. One can 
easily show, however, that the lemma also holds when the elements of 7^12 do not differ in their total fermion number 
by one. We thus conclude that = Hq provided that Na 7^ La ± 1- 

As an example consider the square lattice with doubly periodic boundary conditions as the original graph A. The 
lattice A3 is shown in figure [l] We find that La = 2iVA and the total number of sites in A3 is = 2La + A^a- 
Consequently, this lattice has one ground state at 1/5 filling and one at 2/5 filling. 

Before we continue, a few remarks about the possible boundary conditions are in place. If we have periodic 
boundary conditions on the graph A3 there is no ambiguity, however, there are a few inequivalent possibilities for 
open boundary conditions. On the one hand, we can consider the case that the original graph A has open boundary 
conditions, constructing the A3 graph is then unambiguous. On the other hand, we can have periodic boundary 
conditions on the original graph, construct the A3 graph and only then cut various links to impose open boundary 
conditions. In the A3 graph, however, there are inequivalent links; there are links that connect two S2 sites and links 
that connect an Si site to an S'2 site. In principle we want to allow cutting both types of links, but one has to be 
careful, because it has consequences for the solution to the cohomology problem. When open boundary conditions 
are imposed by cutting links that connect S'2 sites, the cohomology in the sector with La fermions vanishes. This is 
because the S'2 lattice now contains disconnected sites, so the configuration with all Si sites empty is no longer an 
element of 'Hq2 ■ Note that the chain with open boundary conditions and length a multiple of three is an example of 
a A3 graph with this type of boundary conditions. In section [ill A| we indeed found that it has only one zero energy 
state. Conversely, when open boundary conditions are imposed by cutting the other type of links, the configuration 
with all Si sites occupied is no longer an element of 'Hq2 ■ Consequently, for that type of boundary conditions there 
is no ground state with A'a particles. 



Let us now introduce the staggering parameter y in these models as defined in ([2| . First consider the limit of infinite 
staggering. To zeroth order in 1/y the Hamiltonian reads H/y^ fa H2 — {Q2,Q2\- The two elements of Hg, , one 
with Si completely empty and one with Si completely filled, are precisely the ground states of this Hamiltonian. In 



B. Ground states to zeroth order 
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the state with Si completely filled the particles force all S2 sites to be empty. Consequently, this state is annihilated 
by both Q2 and Q2- The state is a simple product state, so the correlation length is finite. We call this charge ordered 
phase a Wigner crystal. 

Now consider the other state which has Si completely empty. The 5*2 sublattice is a collection of two-site chains. 
The ground state of a two-site chain contains one particle that resonates between the two sites: \tpo) = {c\ — c\) / ^/2\0) . 
It follows that in the state with Si completely empty there is one particle on each two-site chain of S2 sites that 
resonates between the two sites. So also this state is a product state. For obvious reasons we call this gapped phase 
a VBS. 

It is interesting to note that although both states are quite simple, they would have been hard to find without the 
cohomology solution. Clearly, the charge ordered state is the expected ground state when the system is perturbed with 
a large conventional staggering that simply favors the Si sites to be occupied. We thus expect this state to dominate a 
significant part of the phase diagram also when the parameters are such that the model is no longer supersymmetric. 
The VBS state, on the other hand, is the energetically favored state when a large conventional staggering favors the 
S2 sites to be occupied. It is remarkable that for the staggering that preserves supersymmetry the two states coexist. 

We now turn to the other limit. For y = we have Q = Qi and we should solve the cohomology problem by looking 
at directly. Since S'l is a collection of disconnected sites, the non-trivial cohomology elements of Hq-^ are all 
configurations which have at least one neighboring site of each S'l site occupied. There are many such configurations, 
with particle numbers in the whole range between N\ and L\. So in the small staggering limit we find that there 
is a large ground state degeneracy. To reveal the nature of the ground states for small but non-zero y we have to 
do degenerate perturbation theory. Remember that the degeneracy is lifted at second order in perturbation theory 
provided that the 'tic-tac-toe' lemma holds. We easily verify that this is indeed the case: since all elements of Hq^ 
have the same number of particles on S'l, namely zero, we must have that Jig — 'Hq^i'HQi)- It follows that the states 
with particle numbers between, but not equal to, Nj^ and La, have an energy of order . We can thus focus on the 
sectors with TVa and particles. In the next sections we construct the effective Hamiltonians for these sector and 
show that the ground states are quantum liquids. 



C. Projected four-state clock model 

In this section we will construct the effective Hamiltonian and obtain an exact expression for the ground state in the 
sector with TVa particles. For simplicity wc first consider the case where the original lattice A is the square lattice. As 



for the chain (see section III B I we obtain the following general expression for the effective Hamiltonian that couples 



the states \i})k) that have zero energy for y ^ Q: 

m 

For F — N\ the states l-^fc) are such that all Si sites have precisely one occupied neighbor. It follows that each site 
can be in four states corresponding to the position of the occupied neighbor. We denote these four states \a) as up, 
down, right and left: \U), \D), \R) and \L) (see figure [2]). With this notation the effective Hamiltonian takes the 
simple form 

^eff,Ar, ^ _ ^ ^P|a)(a'|P+ ^ P^,), (18) 

where P projects onto the states \ipk)- The first term represents flipping the state of the vertex from one orientation 
to another. Due to the hard-core character of the fermions not all flips are allowed, for instance a vertex cannot flip 
into the up-state if the vertex above it is in the down-state. The projection operators P impose these constraints. 
The second term is purely diagonal. The easiest way to compute it is by realising it counts the number of sites on S2 
whose neighbors are empty. We will discuss the second term in more detail when we consider the general case. For 
now we just mention that for the doubly decorated square lattice with doubly periodic boundary conditions this term 
is equal for all {tpk), so it is just an overall constant. 

Remarkably, the first term in the Hamiltonian is precisely the Hamiltonian of the constrained clock model considered 
in [22] apart from an overall constant. In [21] the model arrises as an effective Hamiltonian for ultra cold bosons in 
a strongly tilted optical lattice. The lattice is a singly decorated square lattice which is tilted along the diagonal. In 
the limit of a strong tilt, the true ground state of this system simply has all bosons occupying the site in the lowest 
corner. However, if one starts from a Mott state (one boson per site) and adiabatically tilts the lattice the true ground 
state cannot be reached, because there is too much excess energy in the system. Instead the resonant subspace turns 



11 



□-•-o-n 



□ ^ □ 



□-o-«-n 



□— ► — □ 




FIG. 3: Mapping of elements of the cohomology of Qi at = La to states on the original square lattice, where the 

degrees of freedom are arrows living on the links. 



out to be presicely that of the constrained four-state clock model. It may be surprising that both a bosonic and a 
fermionic system map to the same effective model. This can be understood by noting that in the limit of j/ — > the 
fermions get confined to the sites neighboring one vertex. Consequently, they cannot hop past eachother any more 
and effectively lose their fermionic character. 

The effective Hamiltonian, (18), is again of the Perron- Frobenius type, which ensures that the ground state is 
unique and nodeless. It was shown in 22J that the ground state of the constrained clock model with doubly periodic 
boundary conditions is the equal superposition of all states IV'fe)- This results carries over directly to the ground state 



of ( 18 ) on the doubly decorated square lattice with doubly periodic boundary conditions, since in that case the two 
models only differ by an overall constant. In the supersymmetric model the second term ensures that the ground 
state has zero energy. 

We now turn to the general case and show that, contrary to the constrained clock model, the second term also 
ensures that the equal superposition state is always the ground state independent of the original graph A and the 
imposed boundary conditions To see this we define K2, which counts the number of links between two empty ^2 
sites and ki which counts the number of empty 52 sites at the boundary of the graph. Note that ki is nonzero only if 
we impose open boundary conditions by cutting links in the A3 graph that connect two S2 sites. Using this definition 
we can write the second term in the Hamiltonian as A^a + 2k2 + ni ■ We now proceed as in [2^ , that is, we count the 
number of states to which the Hamiltonian connects some state \ipk)- Note that the first term will change the state of 
the site by occupying a different neighboring S2 site, however, this can only happen if the ^2 site has all neighboring 
sites empty. The number of such S2 sites that can be occupied is precisely 2k2 + Hi ■ Using this result we show that 
each row in the effective Hamiltonian sums to zero, which implies that the equal amplitude state has zero energy. We 
just showed that the off-diagonal elements in each row sum to —2k2 — ki- One easily checks that the diagonal term 
is — A'a + Aa + 2k2 -I- Hi- Clearly, the total sum is zero. Since we know that there is a unique zero energy state for 
y > 0, we find that the equal amplitude state is the ground state. Note that this is also true if the sites in the graph 
A have a different connectivity, this just means that there is a different number of clock states jcr). The fact that this 
holds for general A3 graphs and boundary conditions is again a magical property of supersymmetry. 

In [22] it is shown that correlation lengths in the ground state of the constrained clock model are finite, which 
suggests that the phase is gapped. This result carries over directly to the staggered supersymmetric model, where the 
original graph is the square lattice with doubly periodic boundary conditions. For the other cases, however, it is not 



The only type of boundary wo do not allow is the case where open boundary conditions are imposed by cutting links in the A3 graph 
that connect sites of the different sublattices. As explained before, for this type of boundary condition the cohomology vanishes in the 
sector with N/^ fermions. 
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clear whether the correlation lengths are still finite. Nevertheless, since the ground state is an equal superposition 
state, this can again be investigated using transfer matrices. 



D. Projected two-state link model 

We now turn to the case with F = Ljy. Remember that for y — >■ oo the ground state in this sector is the product 
state of two-site chain ground states. The ground state of a two site chain is the sing let (4 - 4)/^|0). The product 
state thus has singlets resonating on the two ^2 sites that are added to each link of the original graph. We denote 
this state by |VBS) inspired by the obvious similarity with the valence bond solid states. We now show that in the 
limit of 2/ — > the ground state is the state P|VBS), where P again projects onto the states IV'fe)) i-e. the cohomology 
elements of Qi. Before we write the effective Hamiltonian for this sector, we map the configurations \ipk) to arrow 
configurations on the original graph A, where the arrows live on the links. The number of particles in this sector is 
precisely the number of links in A. We replace the two added sites by an arrow pointing towards the vertex which 
has an occupied neighbor on that link. The mapping is shown in figure [3] It follows that all links have an arrow. 
Furthermore, each vertex has at least one arrow pointing inwards. If we now denote the two possible arrow orientations 
on a link by |r), we can write the effective Hamiltonian as 

links r,r' vertices 

where |1) denotes the vertex configuration with 1 arrow pointing in and the rest of the arrows pointing out. As before 
we count the number of configurations to which the Hamiltonian connects some state \4^k)- Note that all arrows can 
be flipped except if that arrow is the only inpointing arrow at some vertex. So the Hamiltonian connects to La — I'l 
states, where is the number of vertices with configuration |1) in the state \'ipk)- One easily checks that the diagonal 
term is precisely this number. To check that the state P|VBS) is a zero energy eigenstate of this Hamiltonian, we 
note that configurations which map into each other by flipping one arrow, have opposite signs in this state. From this 
one readily checks that P|VBS) is a zero energy eigenstate and thus the ground state in the limit of y — 0. 

The projected VBS state is also a quantum liquid state. This can be seen by absorbing the minus sign in the single 
particle singlet states into the deflnition of the creation operator of the fermions. If we define the fermion creation 
operators on each link as fl = c\ and /| = — 4, the single particle state reads (// + /|)/a/2|0) and the projected VBS 
state clearly is an equal amplitude superposition state. We note that the construction of this state is reminiscent of 
the construction of the AKLT state. In the case of the AKLT state, however, one projects out many more states per 
vertex than we do here. 



E. Summary of the different phases 

We conclude that we have found exact ground states in the limits for y — >■ oo and y — for the A3 graphs. In this 
last section we discuss the character of the corresponding phases. In the limit of y — )■ 00 the ground state a.t F = N\ 
is a simple Wigner crystal with all sites of the original graph A occupied, whereas the ground state at F = L\ is the 
VBS state defined above, i.e. singlets on all added sites. It is clear that the correlation length is zero for the state 
with TVa particles and at most two lattice spacings for the state with L\ particles. In the other limit, we obtain two 
types of quantum liquid states. These states are equal superpositions of all states in the low energy subspace. For 
these states we do not know what the correlation length is in general. Both states can be obtained by a projection 
of single particle states onto the space of cohomology elements. This projection induces some correlation. However, 
intuitively one would not expect that this will induce a diverging correlation length in a two dimensional system, since 
it is just not constraining enough. This reasoning suggests that the spectrum in the limit of small y is also gapped in 
both particle number sectors. As we will see in the next section, this is not true in one dimension. 

The ground states in the various gapped phases are quite different in character. In particular, for the case with 
F — Na, where in one limit we have the liquid state and in the other limit the completely ordered state. Nevertheless, 
none of these states breaks any symmetry of the Hamiltonian. This means that there is no obvious order parameter 
we can write down that distinguishes the phases. It follows that it remains unclear whether or not there is a phase 
transition at intermediate values of y. 
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V. MASSLESS KINK IN STAGGERED CHAIN 



In this section wc return to the staggered chain discussed in section [III A[ It is clear that the staggered chain faUs 



into the class of staggered A3 graphs, where the original graph A is also a chain. In section III A we considered open 
boundary conditions imposed by cutting a link that connects two 5*2 sites and found that there is a unique ground 
state. In the limit of large staggering it corresponds to the state with all Si sites occupied. In the limit of small 
staggering the ground state is the quantum liquid state: an equal superposition of all configurations in the low energy 
subspace. In the previous section we argued that for the two dimensional graphs all excitations over the quantum 
liquid state are gapped. Here, however, we will see that the liquid state supports a gapless mode. 

To appreciate the character of the gapless mode we briefly review some known results for the supersymmetric model 
on the chain. For y — 1, the uniform chain, it has been shown that the system is integrable and the continuum limit 
is described by an A/" = (2,2) supcrconformal field theory with central charge c = 1 [l][51[T5]. This indicates that 
the model is gapless for y — 1. Furthermore, it was found that when periodic boundary conditions are imposed and 
the length of the chain is a multiple of three, there are two zero energy states. Staggering the Hamiltonian away 
from the critical point by setting y 1 has been identified with perturbing the continuum theory with the relevant 
supersymmetry preserving operator of conformal dimension 4/3 |161 117] . For the chain with periodic boundary 
conditions and length a multiple of three it was shown that the system flows to a gapped phase both for y > 1 and 
y < 1. Furthermore, it was found that the ground state is twofold degenerate for all values of y, including y = 0. Note 
that this is a special case, because we typically found that the ground state degeneracy of A3 graphs is increased at y 
strictly zero. We understand this by noting that on the periodic chain there are only two configurations that satisfy 
the condition that all sites of the original graph have at least one neighbor occupied: 

• □ o •□ o • • • • □ o •□o, 

oD • oD • • • • o □ • oD • . (20) 

Note that the two ground states are characterized by two types of order. By left-order we denote the phase in which 
the sites to the left of the Si sites are occupied and by right-order we denote the phase in which the sites to the right 
of the Si sites are occupied. 

In this section we find that the gapless mode, supported by the open chain, corresponds to a kink separating the 
two types of order that characterize the two ground states on the periodic chain. Even though this kink is massless, 
the spectrum of the periodic system is gapped because the anti-kink is gapped. This is easy to see in the limit of 
small y; one cannot have the right-order to the left of the left-order in a state where all Si sites have at least one 
neighbor occupied. An anti-kink thus looks as follows 

• • • o □ • oD • oD o •□ o •□ o . . . (21) 

and costs an energy of order 1 in the small y limit. In a periodic system one can only excite kink-anti-kink pairs, 
whereas in an open system one can have one without the other. This explains the somewhat counter intuitive fact 
that changing the boundary conditions changes the system from gapped to gapless. The way one should really think 
about this, is that by changing the boundary conditions one goes to a different sector. In one sector the number of 
kinks equals the number of anti-kinks and in the other sector there is one more kink than there are anti-kinks. 

We will now analyse the staggered chain with open boundary conditions in more detail and show that the kink 
spectrum is similar to the spectrum of a particle in a box. We thus say that this is a one dimensional system which 
supports a zero dimensional massless mode. 

Let us introduce the following notation. If we define 

EE n 4,10) and |i?:r) EE n ci^^m, (22) 

k—n k—n 

we can write the states iV'z) and defined in ([7| and ([s]) as 

i^.) = ii-i)ii?r'), 
!</),) = iL^)ii?r'). 

With this notation the two ground states of the periodic chain of length L = 3j read |i^o~"^) and |i?^~^). It is clear 
that the states |V'i) and are characterized by the location of a kink between the two types of ordering present in 
the two ground states of the periodic chain. 



As we have seen the ground state degeneracy at y = is lifted by the effective Hamiltonian, W (16 1. It can easily 



be verified that this is simply the one dimensional version of the constrained clock model discussed in the previous 
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section. The eigenvalues of W are 



e. = 4sin^(— j, (23) 



with < m < j. The corresponding wave functions arc found to be 



E 

A:=0 



, 7rm(fc + l/2)\ 

cos( ^— Y^llVfc)- (24) 



The ground state (m = 0) is the state in which the kink between the two types of order is completely delocalized. 
Clearly, the kink spectrum is massless, because for large chain lengths L = 3j the energy of the low-lying states 
scales as jj^. Furthermore, we note that the ground state has long range correlations which are associated to 
the delocalization of the kink. To see this we note that in Wigner crystal type ordered states correlators such 
as (non^i) — {no){n3i) would be stricly zero. For example, in the states ipj this correlator vanishes. However, after 
delocalization of the kink due to the gapless mode this correlator becomes non-zero. Indeed, we find for the normalized 



ground state (eqn ( 24 ) with to = 0) 



{nonsi) - (no)(n3i) = 



If we define x = i/j, and remembering that j is proportional to L, we find that this correlator behaves as (1 — x)/L. 
We associate the 1/L decay with the massless kink, which is essentially a zero dimensional mode. It is clear that this 
is in contrast with the exponential decay of correlators observed for the constrained clock model. The fact that we 
do not observe powerlaw decay with x comes from the absence of truly one dimensional critical modes. 

Figure |4] shows a plot of the low lying energies for various chain lengths for y ~ 1/100 alongside the expression in 
eqn ( |23[ ). Very good agreement is observed between the numerically calculated values and the prediction from the 
perturbative expansion. We have also verified an excellent agreement of the wave functions and their perturbative 



expression (24). Note that in the continuum limit both the energies and the wave functions are very similar to those 
of a particle in an infinite potential well. The differences are due to the fact that the kink experiences an attractive 
potential at the edges of the well. We say that the kink spectrum corresponds to a zero dimensional massless mode 
in a one dimensional system. This is in contrast to the massless, quantum critical phase observed for y = 1, which is 
truly one dimensional. 

In the following paragraphs we discuss the density and entanglement properties of the kink phase (see also [29]). 
The numerical results presented here and in the following are from calculations performed using large scale iterative 
exact diagonalisation techniques on large parallel distributed memory machines. The freely available code named 
DoQO (Diagonalisation of Quantum Observables) was employed for this purpose [50] . 



A. The ground state to higher order in y 



Now that we know the ground state to zeroth order we can go back to (11) to obtain the ground state to first 
order in y. As we noted before 7?i only has non-zero elements between configurations which differ in their number of 
particles on the sites 3z -I- 1 by one. Since the jV'fc) have no particles on these sites, i/i maps them to configurations 
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FIG. 4: Plot of the energies of low lying excitations of open staggered SUSY chains with staggering y — 1/100 and 
lengths 3k with 4 < A; < 8. x = 2Ff2 ■ The equation used for the fit is /(n) = 4y^ ^^^'^(wT^^ 



with one particle on this sublattice: 

3 



k=0 



p=0 



= E 

fe=0 

q—k 
j r k—2 

k=0 Lp=0 
q=k+l 

i-1 



+2j2\Lt')®\°'°)k®\Rir+\) 



k=0 



Here we introduced the notation 



o"o)fe=4.+J0). 



(25) 



(26) 



It is easily checked that Hq on all these configurations gives one. It follows that in (11) we have E„i — 1 and we find 

= -Hi|1'o). (27) 

To obtain the wave function to higher order is more subtle, because (^I'ol^'s) will typically not be zero for s > 1. 
This is nicely explained in fTT). The solution is to redefine the ground state as 



\^)=Y,akiy)m+y\^i)+y'\^2) 



(28) 



fc=i 
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with (■0fc|4's) = for all s > and 1 < fc < j + 1. As is explained in [T7], the coefficients ak{y) are polynomials 
in even powers of y. This is directly related to the fact that Hi maps between states with even and odd number of 
particles on the sites 3« + 1. The coefficients of in the polynomials ak{y) follow from the conditions (-0^1 ^'4) = 0. 
Consequently, one has to go to fourth order in perturbation theory to obtain these coefficients. Using numerics we 
compute ak{y) to second order for j up to 7 and, for a convenient normalization, wc find 

ai^aj+i = l + 0{y^), a,n = 1 - + 0{y^) for 2 < m < j . (29) 



B. Density 

The density on site m is given by the ground state expectation value of the number operator (n,„). It is easy to 
see that the density only contains even powers of the staggering parameter. We investigate the density distribution 
in the limit of small y. To zeroth order the ground state is given in 
sites 3z + 1 is zero. Furthermore, one finds that the density on the ot : 
edge. For a chain of length L = 3j it is given by: 



17). Clearly, to zeroth order the density on the 



rer sites depends linearly on the distance to the 



("m) = !!f. , "\ m = 3i 
Hj + 1) 

{nrn) = for m = 3i + 1 

/ \ m+1 
{n„, ) = TTT-. TV tor m = 61 



3(J + 1) 



The density to second order can be obtained from the ground state wave function to first order (25)-(27). Upon 
inspection one finds: 

/ ^ 3j-m , (3j~m)(j-l) + 6(j + l) ^ 
" ^ y 7T7, for rn — 61 



2 



3(j + 1) 3(.? + 1) 

y2 /-iJL for m = ii + l (30) 
(j + 1) 

m + 1 2("'-+l)(j- l) + 60' + l) 



y for TO = 3i + 2. 



^ 3(j + l) ^ 3(j + l)2 

To get the fourth order correction to this result, one has to go to third order for the ground state. This was done 
numerically for chains of lengths up to L = 21. From these results we infer that the density on the sites 3i + 1 to 
fourth order is: 

y^'^lM^ form = l,L- 




^ ^^-y4 3,-+io.,+ii form = 3* + l and l<z<j-2. ^ ' 

Figure [5] shows numerically calculated values of the one point functions along side the values obtained from first 
and fourth order perturbation theory. We note that these results are quite different from the density results obtained 
for the staggered periodic chain [ini [U] . For the periodic chain, the density turns out to have the special property 
that it is scale free. This means that the first 2j terms in the series expansion of the density in y are independent of 
the system size for systems of length L = 3j. Clearly, this special property is absent for the density in the ground 
state of the open chain. This is not surprising since the kink state is very different from the ordered ground states of 
the periodic chain. In the limit of large y, where the ground state of the open chain is very similar to one of the two 
ground states of the periodic chain (both have all Si sites occupied), it is not so clear whether or not the scale free 
property of the density would survive on the open system. We have studied the staggered density numerically in this 
limit for the chain with open boundary conditions and did not find any scale free properties. 



C. Entanglement entropy 

We consider the entanglement entropy in the limit of small staggering parameter y. In this limit the ground state is 
given by the expression in equation |17[ From this we determine the entanglement entropy and entanglement spectrum 
between the subsystems when the chain is partitioned. Any state \^) on the full system can be expressed as a sum 
IV') = J2a p ^a,p\^A)\''l'B) whcrc Wa.p IS Called the weight matrix and \ip\) and IV's) ^.re orthonormal states spanning 
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FIG. 5: Density on staggered sites (n^) for i mod 3 = 1 for (a) y = 1/100 and (b) y ~ 1/10. Shaded plotmarkers 
correspond to the predictions from perturbation theory. In (a) this is to first order using expressions in equation 
(30) and in (b) this is to fourth order using expressions in equation (31). Here we foUow the notation in |31j and use 

fc± = (L ± l)/2 and L = L/3 + 1. 



spaces and H g respectively. When an open chain with 3j sites is divided into subsets A and B containing — 31 
and Ib — 3(j — /) sites, respectively, the weight matrix Wik has m rows and n columns where m^l + l,n = j — l + l. 
The elements of this matrix are: 

Wik = , . , {So,i + Sn-l,k - So,iSn^l.k) 

VJ + 1 

Multiplying this matrix by its transpose results in an m x m matrix which is the reduced density matrix of the subset 
A. 



PA = WW^ = TrB{pAB) 
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Likewise pB = W^W. The matrix elements of pA are: 



^ for i > or fc > 
for i = A; = 



This is a rank two matrix with eigenvalues given by: 

^ _ 1 (m + n-l ^ ^/(m - n - 1)^ + 4(m - 1) ' 

and unnormalised eigenvectors of the form: 



. ^ .n-m + 1 , J{m-n~lY + Aim - 1) 
F±) = I 7, ± 7, 



(32) 



When the chain is divided into subsets with lengths I a = 3Z — 1 and Ib = 3(j — ^) + 1 the weight matrix has 
dimension m x n with m = l + l,n = j — 1 + 2. The elements of this matrix are given by the expression: 

Wik = . . , ^ i5o,i + Sn-i.k - 2(5o,i5„-i,fc) 
V J + 1 

The related reduced density matrix resulting from the multiplication of this weight matrix by its transpose again has 
dimension m x m but in this case has elements given by: 



\PA\ik 



j jij- for i > and fc > 
fori = fc = 



This is again a rank two matrix but with eigenvalues given by: 

^ m — 1 I 71 — 1 J — Z + 1 

""TTT^JTT' j + i 

and unnormalised eigenvectors of the form: 

Ixo) = |1,0,...,0), = |0,1,...,1) (33) 

The entanglement entropy in terms of the eigenvalues of the reduced density matrix is Sy-^ = — log(Ai) and 
the values of the entanglement spectrum are = — log(Ai). Using the expression for the entanglement entropy the 
analytical results were compared to numerical data for an open chain with L = 39 and small value of staggering 
parameter y = ^ (figure [g]) . There is very good agreement between the analytical and numerical data. 

In the limit of infinite system size where the system is partitioned into two equally sized subsets it can be seen from 
the expressions for the eigenvalues that both converge to 1/2. In this case the entanglement entropy is Svn = log(2) 
and the entanglement spectrum consists of two points at log(2). 

The two points of the entanglement spectrum in each case can be interpreted in terms of the kink picture. In the 



case where the subset length is not a multiple of three, it is clear from the eigenvectors (331 that the two eigenvalues 
of the reduced density matrix correspond to the cases where the 'kink' is in subset A only or in subset B only. In 
the case where the chain is split into subsets with lengths that are multiples of three it is possible that as well as 
the kink being in either subset it can also be split by the partitioning such that half is in each subset. We define the 
vector \xa) which corresponds to the kink being completely in subset A and the vector jcc^) which corresponds to the 



kink being either not in A or split by the partitioning. In terms of these vectors the eigenvectors of pA ( 32 ) can be 
expressed as: |a;±) = ^(I^^a) ± l^^^))- In the infinite limit where the eigenvalues are degenerate, the superpositions 

■^{\x^) ± |a;_)) are also eigenvectors. Then + \x-)) — \xa) which corresponds to the kink being completely 

in A and -^{\^+) ~ l^;-)) — \xai which corresponds to the situation where either the kink is not in A or is split by 
the partitioning. 

Finally, we mention that the kink picture of the ground state allows us to derive an expression for the entanglement 
entropy as a function of subsystem size I a by writing a (trivial) wavefunction for the kink in the continuum limit. 
We note that this description does not capture the clear Z3 substructure present in the finite system. Since the kink 
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FIG. 6: Numerically calculated entanglement entropy data for the ground state of the open i = 39 staggered chain 
with staggering parameter y — plotted alongside the analytical results for the small staggering limit indicated by 
the black crosses. The dotted line corresponds to the result for the system in the continuum limit where the Z3 
substructure vanishes. As discussed in the text this line has a very simple interpretation in terms of the kink. 



is distributed uniformly over the chain, the continuum limit wavefunction for the kink is just 5'kink(a;) — 1/v'i for 
< X < L and zero otherwise. Since the kink is either in part A or in part B, the reduced density matrix is 

PA = Pa\xa){xa\ + Pa\xa){xa\- 

Where integrating out part B simply gives the probabilities of finding the kink in part A or part B, which follow 
from the wavefunction above: Pa — Ia/L and P-^ — {Ia ^ L)/L. Plugging this into the formula for the entanglement 
entropy, one finds: 

SvnUa) = -lA/L\og[lA/L] - {L-lA)/L\og[{L-lA)/L\. (34) 

The fact that in this case the levels of the entanglement spectrum can be so directly related to a physical feature of 
the system, namely the kink is a very interesting feature of this model. 



VI. OUTLOOK 



An obvious follow-up on this work is to study the staggered supersymmetric model on other lattices, in particular 
the ones that exhibit its key feature: superfrustration. In [2] the limits of small and large staggering were considered 
for the supersymmetric model on the martini lattice. The martini lattice has an exponential number of ground states 
and the ground states are in one-to-one correspondence with dimer coverings of the underlying honeycomb lattice. 
In one limit of the staggering parameter the dimer coverings become exact ground states. It was shown that by 
introducing another parameter that breaks supersymmetry, this ground state degeneracy is lifted and the effective 
Hamiltonian is a quantum dimer model. 

We will briefly discuss two systems that seem particularly promising for continuing the strategy introduced in this 
paper in which supersymmetry is preserved: the square lattice decorated by one extra site on each link and the 
octagon-square lattice, which is obtained from the square lattice by replacing each vertex by a smaller square that 
is tilted by 45 degrees. When wrapped around a torus, both lattices show an exponential ground state degeneracy. 
Still the frustration is suppressed, since the ground state entropy is subextensive: it grows with the linear size of the 
system, not the area. 

Let us first consider the octagon-square lattice. The cohomology problem was solved in [2]. Interestingly, the 
cohomology problem can be solved for various choices of the two sublattices and 82- Clearly, this makes this an 
interesting system to look at, since various staggerings can be considered. Let us briefly discuss a choice that seems 
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particularly promising. Let the top site of each added square be in Si, the S2 is a collection of chains with length 
a multiple of three. If periodic boundary conditions are imposed in the horizontal direction, that is in the direction 
along the chains, the elements of Hq^ have all 5*1 sites empty. Since each periodic chain has two ground states, the 
number of elements is 2^" , where Nc is the number of chains. On the other hand, the elements of Hq^ also have all 
Si sites empty, however, now the particles on the ^2 sites, must be arranged such that all Si sites have an occupied 
neighbor. It follows that all chains can be in one of the two configurations depicted in (20). So again the number of 
elements in 2^"=. So even at strictly zero staggering, the number of zero energy states is unchanged. It follows that 
one cannot write an effective Hamiltonian that couples these states. What is interesting, however, is that for very 
large and very small staggering, this system behaves as a collection of one dimensional systems, which in one limit are 
quantum critical and in the other limit are gapped. If you wish, there is some resemblance with the one dimensional 
chain, where we had a zero dimensional gapless mode in a ID system, whereas here we have one dimensional gapless 
modes in a two dimensional system. The question is of course what happens to the criticality when the staggering 
parameter is finite. Since the critical chains are sensitive to a relevant perturbation that drives them to the gapped 
phase that we observe in the limit of very small staggering, a possibility is that the system is gapped for all finite values 
of the staggering parameter. In that case the question arises as to whether these 2^" ground states are topologically 
distinct or not. 

Another interesting system to look at is the square lattice with one extra site on each link, which we will call the 
singly decorated square lattice. This lattice falls into the class of A2 graphs. As we mentioned it can be proven 
[m [28] that the cohomology problem of these graphs maps to the cohomology problem of the original graph via a 
particle-hole inversion. The cohomology problem of the square lattice with doubly periodic boundary conditions has 
been solved [TT]. The ground states are found to have densities v = F/N, where TV is the total number of sites in the 
whole range between 1/5 and 1/4. For the singly decorated square lattice this means that there are zero energy states 
for all rational densities in the range [1/4, 4/15]. This is true for all non-zero values of the staggering parameter. 
As before we choose the sublattice to be the original lattice and 5*2 to be the added sites. For the staggering 
parameter strictly equal to zero, all cohomology elements of Hq^ are zero energy states. It turns out that these 
states can be found for all densities in the range [1/6, 2/3]. At 1/6 density, the cohomology elements are one-to-one 
with close packed dimer configurations on the links of the original square lattice. For any non-zero staggering we 
know that there are no zero energy states at 1/6 filling. Interestingly, we find that to second order the effective 
Hamiltonian in this sector is just a constant, so the states gain some finite energy as they should, but they remain 
degenerate. This degeneracy is lifted at fourth order, where the effective Hamiltonian is the usual plaquette flip term 
in the Rokhsar-Kivelson model [H]. Since the potential term is absent (it arises only at eighth order) the system is 
in the gapped phase which shows columnar ordering of the dimers. The study of the effective Hamiltonians in the 
sectors at higher densities is beyond the scope of this project. From the cohomology structure is it clear that they 
may be quite involved, but therefore possibly also quite interesting. 



VII. CONCLUSION 



We showed that exact ground states for the supersymmetric model on decorated graphs can be obtained in the 
limit of small and large staggering. In one limit we find that the ground states are either a Wigner crystal or a simple 
VBS, in the other limit the ground states are quantum liquid states arising from a projection of a product state with 
one particle per plaquette (for F = Na) or link (for F = La), reminiscent of the AKLT quantum liquid states. It 
turns out that the effective Hamiltonian in the sector with F = N/^ for the decorated square lattice with periodic 
boundary conditions is precisely the constrained clock model of [22] . Since the ground state is an equal superposition 
of classical configurations techniques from classical statistical mechanics could be used to study correlation functions 
in this state. They found an exponential decay for all correlation functions considered. This suggests that the system 
is gapped. Intuitively, we expect this to hold for all liquid states that arise on the A3 lattices. The liquid state on the 
chain forms an exception since it does show long range correlations. The massless mode in this system is understood 
as a kink separating two types of order. Apart from studying correlation functions, it would be interesting to look 
at the entanglement entropy of the liquid states, since a saturating entanglement entropy in a finite range of the 
subsystem size is a clear signature of a gap. 

The question of whether or not there is a phase transition for intermediate values of the staggering parameter can 
be addressed by numerical studies of ladder systems. The chain is clearly a very special case, because of the twofold 
degeneracy of the ground state when periodic boundary conditions arc imposed. Decorated ladders, however, are 
expected to behave very similar to full two dimensional systems. 

Finally, a promising follow-up on this work would be to look at other lattices for which the cohomology problem 
can be solved in two limits of the staggering parameter. As we mentioned, the supersymmetric model typically has a 
vast ground state degeneracy, which will make this analysis more challenging. However, it also holds the promise of 
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discovering exotic phases with frustration. 
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